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I. INTRODUCTION 

During the last decade, the relationship between quan- 
tum phase transitions and entanglement has become an 
important research domain^ . Although it is natural to 
expect some deep changes in the ground state of a system 
at a transition point, the real problem is to measure these 
variations or, in other words, to characterize the quan- 
tum state structure. In most cases, the study of an order 
parameter or the behavior of correlation functions is suf- 
ficient to detect and analyze a phase transition but one 
may wonder whether more "intrinsic" measures could be 
helpful. Following pioneering works in one-dimensional 
spin models ^ , many studies have been devoted to this 
problem (see Ref. for a review), but only a few allow 
for an exact solution in the thermodynamic limit which 
is a key ingredient to characterize a phase transition. 

The goal of this paper is to propose a class of sim- 
ple, fully connected (collective) models in which ground- 
state entanglement properties can be studied in details. 
These models, in which degrees of freedom (spins 1/2) 
mutually interact, can be seen as generalizations of the 
Lipkin-Meshkov-Glick model'^'"' for which most entangle- 
ment features are now well known*^^^'\ The main reason 
for introducing these collective systems is that they not 
only allow us to study a second-order phase transition as 
in the Lipkin-Meshkov-Glick model, but also allow the 
study of first-order transitions. As surprising as it may 
seem, although "collective" may be thought to lead to a 
pure mean-field behavior, we will see that the is deeply 
entangled. Furthermore, although it might be naively ex- 
pected that entanglement measures will simply display 
jumps at first-order transitions, we will show that this 
does not hold for one of the models, the spectral prop- 
erties of which show some similarities with those of a 
system exhibiting a second-order transition. 

The study of these collective systems is also moti- 
vated by the fact that they are much simpler to ana- 
lyze than their nearest-neighbor counterparts on a finite- 
dimensional lattice. Indeed, most entanglement measures 



rely on a multi-partition of the microscopic degrees of 
freedom. For instance, the concurrence is obtained by 
separating a system of N spins into two parts (of sizes 
N — 2 and 2), the Renyi entropy is obtained by split- 
ting it into two parts of arbitrary sizes {N — L and L), 
a tri-partition is required to compute the negativity' of 
a mixed state, etc. Thus, the main problem often con- 
sists of computing reduced density matrices for a given 
partition. In the models studied below, this crucial step 
can be achieved since the original spin problem can be 
mapped onto a quadratic bosonic Hamiltonian. 

The structure of this paper is the following. In Sec. II, 
we introduce a family of models and compute their low- 
energy spectrum (ground-state energy and gap). This 
allows us to determine their phase diagram and char- 
acterize the quantum phase transitions. Analytical ex- 
pressions are obtained in the thermodynamical limit and 
compared with exact diagonalization results. In Sec. Ill, 
we discuss the ground-state entanglement by focusing on 
three different measures : the concurrence, the Renyi en- 
tropy, and the negativity, which rely on a one-mode, two- 
mode, and three-mode description of the bosonic Hamil- 
tonian, respectively. Once again, exact results in the 
thermodynamical limit are compared to numerical data 
for a representative set of parameters. 



II. MODELS AND QUANTUM PHASE 
TRANSITIONS 

A. Hamiltonians 

We consider a system made of N spins 1/2 whose 
Hamiltonian reads 



H = -N 



cos w 



, sm uj 



(1) 



In the above equation, Sa — '^jLi'^'j /'^ are total spin 
operators along the a — x, y, z direction, with a" being 
the usual Pauli matrix at site j and S = N/2 denotes 
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the maximum spin value. This class of models is defined 
by the non- negative integer parameters (m, n). In what 
follows, we shall refer to a model with given values of m 
and n as the (to, n) model. 

Without loss of generality, we shall restrict ourselves 
to TO ^ n ^ 1. We furthermore exclude the trivial case 
m = n — 1^ since it describes a large spin in a mag- 
netic field, and only displays a crossover, but no quan- 
tum phase transition. The (2, 1) model is a collective 
version of the transverse-field Ising model, known as the 
Lipkin-Meshkov-Glick model " . The (to > 2,1) models 
are multispin generalizations of such a model. The (2, 2) 
model can be seen as a collective version of the quan- 
tum compass model"''. Note that, in two dimensions, the 
latter is dual to the Xu-Moore model"'""'', but the col- 
lective version of the Xu-Moore model, namely, the (4, 1) 
model, is not dual to the (2, 2) model (in fact, as will be 
seen below, the latter two models have rather different 
properties). 

As we shall see, all models under consideration exhibit 
a quantum phase transition when the control parameter 
uj G [0,7r/2] is varied. This quantum phase transition 
occurs at uj = 7r/4, provided one imposes K^^n to take 
the following value'^'': 

K2,i = 2, (2) 

TO™/2(to-2)™/2-1 

Km>2,l - 7 , N„, 1 , (3) 

(to — 1)™ ^ 

^m^2,n>2 = 1- (4) 

These values can be found easily (see Sec. II B 2 for de- 
tails) . 

Finally, let us note that all Hamiltonians preserve the 
magnitude of the total spin, i. e., [H, S'^] — 0. When to 
(n) is even, the Hamiltonian furthermore has a spin-fiip 
symmetry since it commutes with J^^- cr| (J^^- crj). 



tra between these two figures, one can infer the behavior 
of the various models in the thermodynamical limit. 

The (2, 1) model displays a collapse of levels onto the 
ground state, at the transition point uj = 7r/4, but the 
ground-state energy displays no cusp. From these fea- 
tures, the quantum phase transition is likely to be of 
second order. All other models (except the (2, 2) model) 
have avoided level crossings, which tend to become true 
level crossings in the thermodynamical limit. In particu- 
lar, the ground-state energy has a cusp, and the transi- 
tion is of first order. Note that, in Fig. 2, the darker 
regions are those where one finds many levels (which 
are finite-size precursors of singularities in the density 
of states in the thermodynamical limit). But, contrary 
to what happens for the (2, 1) model "'' ' ', these regions 
do not touch the ground-state energy at the transition. 
This is, however, not true for the (2, 2) model, which dis- 
plays both a collapse of levels onto the ground state as 
well as a cusp in the ground-state energy. From the sec- 
ond feature, one concludes that the transition is first or- 
der, although the first feature is reminiscent of a second- 
order quantum phase transition. Let us stress that 
the (2, 2) model is trivially integrable at the transition 
point, since its Hamiltonian reads H = — ^^(S'^ — Sy), 
where S"^ = {N/2){N/2 + 1). This additional symme- 
try {[H,Sy] = 0) is responsible for the presence of non- 
avoided level crossings at the transition point, even at 
finite N. Actually, all these models are exactly solvable, 
but not in such a trivial way '" ' '. 

It is this variety of behaviors displayed by the different 
models that motivates the study of entanglement mea- 
sures and their sensitivity to the characteristics of the 
quantum phase transition. Before turning to this, we 
shall, however, provide analytical results for the low- 
energy spectrum, which will allow us to introduce the 
basic techniques needed to perform analytical computa- 
tions of entanglement measures. 



B. Quantum phase transitions 

1. Numerical spectra 

Physically, the quantum phase transition stems from 
the competition between the ferromagnetic TO-spin inter- 
action in the x direction and the ferromagnetic n-spin 
interaction in the z direction if n > 1, or magnetic field 
in the z direction if n = 1. 

A numerical study provides an idea about the phase 
transitions of the various models. Such a study can be 
performed for rather large number of spins since, as al- 
ready mentioned, the Hamiltonians commute with S^. 
The collective and ferromagnetic nature of the interac- 
tions implies that one can focus on the maximum spin 
sector S = N/2 (of dimension -I- 1) where the ground 
state is found. The energies per spin e of six different 
models are shown in Figs. 1 and 2, for A'^ = 16 and 
A^ = 256 respectively. From the evolution of the full spec- 



2. Ground-state energy 

Since a large spin behaves classically, a classical anal- 
ysis is expected to provide exact results for the ground- 
state energy of the collective models in the thermody- 
namical limit. We are therefore led to substitute the 
spin operators by their expectation values, namely 

{{So;} ,{Sy) ,{Sz)) = Y (sin e COS (f>, sin sin 0, cos 9) , 

(5) 

where 6 G [0,7r] and (f> G [0, 27r[ are the usual angles 
of spherical coordinates. They are the variational pa- 
rameters that will be tuned to minimize the associated 
classical energy per spin 

e(0,(/)) = — cos aj(sin6' cos (/))™ — „ sina;(cos (?)". (6) 

When TO > n = 1, the analysis follows the mean-field 
calculation of Ref. '() [see also Refs. V.-],'M> for the (2, 1) 
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FIG. 1: (Color online) Spectra (energies per spin) of six mod- 
els as a function of the control parameter cj, for a system of 
A'^ = 16 spins, in the maximum spin sector 5* — N/2. Left 
(right): n — 1 (n = 2). From top to bottom ; m = 2, m = 3, 
and m = 4. 



FIG. 2: (Color online) Spectra (energies per spin) of six mod- 
els as a function of the control parameter uj, for a system of 
A'' = 256 spins, in the maximum spin sector S = N/2. Left 
(right): n = 1 (n = 2). From top to bottom : m = 2, m = 3, 
and m = 4. 



model]. At "large" w (close to 7r/2), the state (6*0 = 0, 0) 
is the only minimum, with energy e(0, — —Km.i sincj. 
At "small" w (close to 0), 0o = = or tt) when 
m is odd (even) , and the angle 9q minimizing the energy 
satisfies mcos6'o sin™"^ 0o = -S^m.itanoj. Requiring that 
the transition take place at = 7r/4, one is led to solve 
the following system of equations: 

Km,i = mcoseosm'^-'^ 9*0, (7) 
Km,i = sin™0*-f i^^aCOS^S, (8) 

where the second equation stems from the continuity of 
e at the transition, and 9^ is the value of 9o at the tran- 
sition, in the small-oj phase. The solution of this system 
yields Eqs. (2) and (3), as well as cos = It is 

therefore clear that, except for the (2, 1) model, 9o is dis- 
continuous at the transition, which is thus of first order 
for all (m > 2, 1) models. 



For the (2, 1) model, the transition is of second order 
(see e. g. Ref 3G), as can be seen from the discontinuity 
of the second derivative which jumps from the value 
~3V2 at w = {ir/iy to at w = (7r/4)+. The large-w 
phase is a symmetric phase with a non-degenerate ground 
state, while the small-w phase is a broken phase with a 
doubly-degenerate ground state, the broken symmetry 
being the parity Sx ^ —Sx. The validity of the classical 
analysis can be assessed in Fig. 3 where numerical data 
can be seen to converge to the classical result. 

When m ^ n ^ 2, one can proceed in the same way. 
One finds that the transitions are all of first-order nature, 
that Eq. (4) has to hold in order to have a transition at 
Lo = 7r/4, and that the angles 9o and 4>o take the following 
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FIG. 3: (Color online) Ground-state energy per spin eo of six 
models as a function of uo, for N = 16, 32, 64, 128, 256 and 
in the thermodynamical limit (red thick line with dots). Left 
(right): n = 1 (n = 2). From top to bottom : m = 2, m = 3, 
and m = 4. 



values 



ijj 7r/4 


uj ^ 7r/4 


00 = (tt) 


00 = (tt) 
any (j)Q 



(9) 



The values in parentheses are other possible values de- 
pending on the parity of m and n. For w ^ 7r/4, the 
states (7r/2,0) and (7r/2,7r) are degenerate when m is 
even, whereas for a; > 7r/4, the states {Q,(f)o) and (tt, (/)o) 
are degenerate when n is even. The degeneracies are al- 
ready predictable from Fig. 1, and the validity of the 
classical results can again be checked in Fig. 3. 



3. Gap 

In order to conclude the analysis of the spectrum and 
of the quantum phase transition of these models, let us 
now turn to the computation of the gap, in the maximum 




FIG. 4: (Color online) Gaps to the first two excited states of 
six models as a function of w, for A'^ = 16, 32, 64, 128, 256 
and in the thermodynamical limit (red thick line with dots). 
Left (right): n = 1 (n = 2). From top to bottom : m = 2, 
m = 3, and m = A. 



spin sector S — N/2. We follow the procedure described 
in Refs. ' ' and refer the reader to these references for 
details. As a first step, we perform a rotation around 
the j/-axis, in order to bring the z-axis along the classical 
magnetization direction 




^ cos 00 sin 00 \ (sj" 

I \ Sy 

sin 00 cos0o/ 



(10) 



Next, we make use of the bosonic Holstein-PrimakofF 
representation of the rotated spin operators'^' 



~ N . ~ 

S^^ a^a and S+ = {N 



a' a a — 



(11) 



where S± = S^iiSy and the a operator is a bosonic anni- 
hilation operator, satisfying [a, a^] = 1. As we shall only 
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focus on the thermodynamical hmit, it will be sufficient 
to keep terms of order and iV" in the Hamiltonian, 
neglecting all terms that go to zero as -> oo, (assum- 
ing a finite number a^a of bosons). For example, in the 
large-w phase, one can write Sx = Sx = N^^'^{a^ + a)/2. 
Note that no term of order y/N appears, thanks to the 
rotation wc have performed. The Hamiltonian then reads 

H ^ Nen + -i + &a^a + -i {a}"^ + , (12) 

where cq, 7 and 5 have the following expressions (the 
value of 9q has been given in Sec. HB 2) 

Co = — cosw(sin0o)™ ~ ^in,n sinw(cos0o)", (13) 

7 = [m(m - 1) cosw(sin0o)'"^^(cos0o)^ 

+n{n~ I)ii:m,„sina;(cos6'o)""^(sin0o)^], (14) 

5 = 2 [m cos a;(sin ^o)"* + ?^^m,n sin a;(cos0o)"] 

+27. (15) 

Note that eo is simply the minimum of the classical 
ground-state energy (6). 

Such a quadratic Hamiltonian is diagonalized via a Bo- 
goliubov transformation 

a = cosh(e/2)fo + sinh(e/2)6^ (16) 

where 6 is a bosonic annihilation operator, satisfying 
[6,6^] =1. The value of diagonalizing the Hamiltonian 
satisfies tanhG = e = —2^/5. With these notations, 

F = iVeo + 7+^ (Vl-e^-l) +A6t5, (17) 

where the gap is A = 5^/1 — e^. Let us emphasize that 
A is the gap above the possibly-degenerate ground state, 
but does not capture the energy splitting between the 
ground states if they are degenerate. We compare the 
spectrum of Eq. (17) with numerics in Fig. 4 for the gap 
to the first and second excited states, so we get at least 
one (and possibly two) nonzero value in the thermody- 
namical limit. The relevance of this simple "spin- wave" - 
like approach can be appreciated. 



representation. The necessary steps for the computa- 
tion of the concurrence have already been performed in 
Sec. HBS, but let us give the key ingredients that are 
useful to obtain the other entanglement measures. 

The very first step is to perform the rotation (10). 
Then, one splits the system into p subsystems, so that 

the the spin operators read as Sa — Yl^i=i^a\ where 
a = y ov z. Depending on the entanglement measure 
one wishes to compute, one has p = 1, 2, or 3. One 
then introduces p bosonic operators and their conju- 
gates a\ for each subsystem. Denoting the number of 
spins of each subsystem by iV^, with 'Y^\=i = ^ ^ the 
p Holstein-Primakoff representations read as 

5« = ^ - a\a, and = (n, - a\a>j a,. (18) 

One can then insert these expressions in the Hamiltonian, 
expand all operators and keep terms of order A^^ and , 
neglecting contributions that vanish in the thermody- 
namical limit. After simple algebra, one gets a quadratic 
Hamiltonian 

p p 

1=1 k.l=l 

(19) 

where Bq, 7, and S are given in Eqs. (13)-(15) and where 
Ti = Ni/N. Let us note that, to obtain this pre- 
cise quadratic form, with only diagonal boson-conserving 
terms, one must get rid of terms of the form a^o; 
with k ^ I. To this end, one should use the rela- 
tion 5*^ = S{S + 1), written in the bosonic language, 
namely, 

V^(4«i + H-C-) - 2 ^(1 - n)a\a,. (20) 

k^l i 

Of course, Eq. (19) yields Eq. (12) when only one bosonic 
mode is considered. 

In the three subsections that follow, we shall give a 
minimal amount of computational details, knowing that 
these can already be found in the literature. 



III. ENTANGLEMENT MEASURES 
A. Technical prerequisite 

We shall now compute three entanglement measures, 
namely, the concurrence, the entanglement entropy, and 
the logarithmic negativity. These measures have al- 
ready been computed for the (2,1) model in Refs. 11- 
13,15,17,23. As can be inferred from these works, the 
analytical computations require to write the spin opera- 
tors as the sum of one, two or three spin operators, for 
the concurrence, entanglement entropy and negativity re- 
spectively. Then, one should use the Holstein-Primakoff 



B. Concurrence 

The concurrence C measures the entanglement be- 
tween two spins half, these spins being in either a pure or 
a mixed state . Here, we are interested in quantifying 
the entanglement between any two spins, the others being 
traced over. Finding the concurrence amounts to com- 
puting the entries of the reduced density matrix, which 
can be done easily for symmetric states . However, ex- 
cept in the case of systems possessing a spin-flip symme- 
try, finding a simple analytical formula for the concur- 
rence is not such an easy task'''. Although we have no 
formal proof, we have checked numerically for finite-size 
systems and a couple of values of m and n (even when 
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FIG. 5: (Color online) Concurrence of six models as a function 
of w, for A'' — 16, 32, 64, 128, 256 and in the thermodynamical 
limit (red thick line with dots). Left (right): n = 1 (n = 2). 
From top to bottom : m — 2, m = 3, and m — 4. 



m and n are odd so that there is no spin-flip symme- 
try), that the rescaled concurrence Cr could be simply 
expressed as 



(A^-1)C==1 



4(gg) 
N 



(21) 



This rescaling is needed here because each spin shares en- 
tanglement with its iV — 1 "neighbors" . Thanks to the 
Holstein-PrimakofF representation (11) and to the Bogoli- 
ubov diagonalization of the associated Hamiltonian (16), 
one can show that in the thermodynamical limit 



lim 

N^oo 



N 



1-e 



1 



(22) 



where £ is given just before Eq. (17). Figure 5 displays 
numerical results for increasing system sizes which clearly 
converge toward the expression computed above in the 
thermodynamical limit. 

As can be inferred from this figure, the concurrence 
is cusped but continuous at the second-order quantum 



phase transition for the (2, 1) model, while it displays a 
jump at the first-order transition of the (m, n) models, 
except for the (2, 2) model where it shows a cusp and 
is continuous. The spectral pecuharities of the latter, 
which have been discussed in Sec. II B 1, do not lead to 
a discontinuous concurrence. It therefore seems, in this 
very special case, that an entanglement measure such as 
the concurrence is more sensitive to the "level collapse" 
on the ground state (Anderson's tower structure), than 
to the level crossing, when both effects are present. We 
shall show that this conclusion remains valid for the other 
entanglement measures we have calculated, starting with 
the entanglement entropy. 



C. Entanglement entropy 

The ground-state entanglement between two comple- 
mentary subsystems A and B can be quantified by the 
Renyi entropy, defined by 



1 



1 



In [Tr(p^'' 



(23) 



In the above equation, = Trgp is the reduced den- 
sity matrix of subsystem A (p is the ground-state density 
matrix) and g is a positive number. In the limit g — >■ 1, 
one recovers the usual von Neumann entropy, namely 
£ = Imiq^iSg = — Tr[p_4 lnp_4]. The technique for com- 
puting p_4 has been exposed in Refs. 15,1 . To use this 
method one simply needs the Bogoliubov transformation 
which diagonalizes the Hamiltonian (19) for p = 2 given 
in Appendix A. In summary, for subsystems A and B of 
sizes N_4^ — tN and iVg = (1 — r)iV, one has (in the 
thermodynamical limit and in an appropriate basis) 



PA 



cxp 



In ( ^^—^ — ] c^c 



with (24) 
(25) 



^Ji = v/[T + (l-T)/a][(l-T)+r/a], 



where c and are bosonic annihilation and creation op- 
erators and where a has been defined in Eq. (22). It is 
then straightforward to compute the Renyi entropy 



^^[q\n2-\n[{p + lY-{^^-lY]], (26) 



as well as the von Neumann entropy 



(27) 



We have computed the latter numerically. As can be 
seen in Fig. 6, when the system size grows, the numerical 
results converge to the analytical expressions obtained 
above (to which one must in fact add a term In 2 when 
the ground state is two- fold degenerate). 

One can furthermore see that similar conclusions to 
those for the concurrence can be drawn here. Indeed, 
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FIG. 6: (Color online) Von Neumann entanglement entropy 
of six models as a function of cj, for A*' — 16, 32, 64, 128, 256 
and in the thermodynamical limit (red thick line with dots). 
In all cases, the system is separated in two parts of equal sizes 
iVi = iV2 = N/2. Left (right): n = 1 (n = 2). From top to 
bottom : m — 2, m = 3, and m — 4. 



or in a pure state), was already worked out for the (2, 1) 
model" ' and is obtained as follows. The system is divided 
into three subsystems A, B and C of respective sizes Ni , 
N2 and N3. One then traces the ground-state density 
matrix over one of the subsystems, say B, to obtain the 
reduced density matrix = Tr^p. The logarithmic 
negativity C is defined as 



£ = In Tr 



(28) 



where T^ denotes the partial transposition with re- 
spect to subsystem A. In a basis of states !</),'!/') = 
\'1>)a ® W)c^ this operation reads as ((/)', |0, V') = 

(0, '!/''lft4cl0'i V')- Since £ measures the entanglement be- 
tween subsystems A and C, we would obtain the same 
result by considering Tc in Eq. (28). 

In Ref. 23, Wichterich et al. have shown that once 
the Hamiltonian is written as a three-boson Hamiltonian, 
that is (19) withp = 3, the logarithmic negativity is given 

by 



£ = - ^ In 
2 



1+9- v/52+4Tir3(a + l/a-2)l , (29) 



with 



5 - [ti + T-3 - (n - T3)2] {a + 1/a - 2)/2, (30) 

where a is given in Eq. (22). This result, which is valid 
in the thermodynamical limit, is plotted in Fig. 7. One 
can furthermore see in this figure that the finite-size data 
from exact diagonalizations converge to the value (30) 
when the system size N grows. In addition, all that was 
said for the behavior of the concurrence of the various 
models under investigation holds again here for the log- 
arithmic negativity. 



IV. CONCLUSION 



the von Neumann entropy of the (2, 2) model diverges 
at the first-order transition point, like the entropy of 
the (2, 1) model but contrary to the entropy of all other 
models which is finite but discontinuous at the transi- 
tion. In fact, the entropies of the (2, 2) model and of the 
(2, 1) model diverge logarithmically at the transition, as 
(l/2)lniV and (l/6)lniV, respectively (see Refs. 12, ir,). 
So, once again, from an entanglement perspective, the 
peculiar first-order transition of the (2, 2) model looks 
like a second-order transition. 



D. Logarithmic negativity 

As a final study of the ground-state entanglement 
properties of our class of models, let us compute the loga- 
rithmic negativity- ' . This quantity, which quantifies the 
entanglement between any two subsystems (in a mixed 



The concurrence, the entanglement entropy, and the 
logarithmic negativity, although different entanglement 
measures, show similar features when used to character- 
ize the quantum phase transitions of the class of collective 
models we have introduced in this paper. However, when 
the transition is of first-order nature but accompanied by 
a collapse of levels on the ground state [see model (2, 2)], 
as is usually characteristic of second-order transitions, 
the entanglement of the ground state does not show any 
discontinuity at the transition, but behaves exactly as in 
a usual second-order transition. In such a situation, one 
may wonder whether other "intrinsic measures" would be 
more sensitive to this discontinuous transition. One may 
think about studying the fidelity that has already been 
analyzed for the (2,1) (Lipkin-Meshkov-Glick) model at 
zero"-''''^""''^ and at finite temperature ' or to the ge- 
ometric entanglement, which is also known for the (2,1) 
case . We have computed these quantities for the (2,2) 



8 




FIG. 7: (Color online) Logarithmic negativity of six mod- 
els as a function of u, for A'' = 16, 32, 64, 128, 256 and in 
the thermodynamical limit (red thick line with dots). In all 
cases, one first traces over N2 ~ N/2 spins, and the logarith- 
mic negativity is computed for the remaining spins which are 
partitioned in two subsystems of equal sizes A^i — N3 = N/4. 
Left (right): n = 1 (n = 2). From top to bottom : m = 2, 
m = 3, and m — 4. 



model and found that their behavior is similar in the 
(2,1) and the (2,2) cases although, as already underlined, 
finite-size scalings are different. We wish to underline 
that it this not an isolated case since, for all (m, m) with 
TO ^ 2, one has a first-order transition and the symmetry 
of the Hamiltonian under the exchange (x o z) implies 
that entanglement measures must be continuous. How- 
ever, as can be checked from the exact formulas given in 
this paper, there is no entanglement for (to, n) models 
when TO ^ n ^ 3. Thus, from this perpspective the (2, 2) 
model is a bit singular. 

To conclude, let us emphasize that we focused here on 
the ground-state entanglement. Nevertheless, it would 
be worth considering the full spectrum of these mod- 
els to investigate finite-temperature entanglement, which 
may unveil interesting properties ' . This is beyond the 
scope of this paper but it will be the topic of a forthcom- 
ing publication^''. 



Appendix A: Diagonalization of the two-mode 
Hamiltonian (19) 

To compute the entanglement entropy, one needs to 
diagonalize the Hamiltonian (19) for p = 2. This is done 
by performing the following Bogoliubov transformation : 

ai = [cosh(e/2)6i +sinh(e/2)6|] 0^ + 62%/^, (Al) 
02 = [cosh(e/2)6i +sinh(e/2)6|] 62^^, (A2) 

where tanhG = s = —2^/S. New bosonic operators mu- 



bi,b[ 



b-i,b 



1. 



tually commutes and satisfy 

Inserting these relations in Eq. (19) for p = 2, one 
gets : 



(A3) 

where Ai = 5^/1 — and A2 = (5. 
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